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Molecular building blocks interacting at the nanoscale organize spontaneously into stable mono- 
layers that display intriguing long-range ordering motifs on the surface of atomic substrates. The 
patterning process, if appropriately controlled, represents a viable route to manufacture practical 
nanodevices. With this goal in mind, we seek to capture the salient features of the self-assembly 
process by means of an interaction- site model. The geometry of the building blocks, the symmetry 
of the underlying substrate, and the strength and range of interactions encode the self-assembly pro- 
cess. By means of Monte Carlo simulations, we have predicted an ample variety of ordering motifs 
which nicely reproduce the experimental results. Here, we explore in detail the phase behavior of 
the system in terms of the temperature and the lattice constant of the underlying substrate. Our 
method is suitable to investigate the stability of the emergent patterns as well as to identify the 
nature of the melting transition monitoring appropriate order parameters. 

PACS numbers: 89.75Fb, 81.16.Dn, 81.16.Fg, 05.65. +b, 89.75.Kd 
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I. INTRODUCTION 



Molecular building blocks interacting at the nanoscale 
organize spontaneously into stable ordered monolay- 
ers p]|. The patterns that emerge upon self-organization 
of simple supramolecular units display a variety of sym- 
metries and local ordering motifs with different degrees 
of packing on the surface. Currently such surface coat- 
ings are under active research due to their potential ap- 
plications in nanoscience, such as surface functionaliza- 
tion 0, sensor surfaces Q, or molecular electronics Q. 
The theoretical challenge is to provide tools that allow 
to predict the patterns without performing the actual 
measurements for the supramolecular units Q. Most 
approaches rely on atomistic modeling where the stable 
conformations are calculated using experimental meas- 
urements of the molecular conformations on the surface 
as input. There, the molecular subunit is resolved in 
all chemical details, resulting in a large number of de- 
grees of freedom that have to be included in the mod- 
eling. In molecular mechanics calculations (MM) Q 
and molecular dynamics simulations (MD) 0, Hfl the mo- 
lecular interactions are parametrized by semi-empirical 
force fields which requires a significant number of addi- 
tional parameters. In MM the starting conformations are 
energy-minimized resulting in a stable configuration, yet 
the method intrinsically provides only the closest local 
minimum. In contrast, MD simulations sample differ- 
ent minima and incorporate, in principle, also dynamical 
reconfigurational processes. The quantum aspects such 
as the electronic density of states (DOS) measurable in 
scanning tunneling microscopy (STM) can be computed 



using density functional theory (DFT) 0, yet these 
calculations are limited to small systems. 

Recently, a complementary approach has been intro- 
duced employing coarse-grained models which aim at pre- 
dicting certain features of the patterns with the bene- 
fit of a great reduction of complexity. For example, us- 
ing an effective hamiltonian accounting for the energetics 
of the respective orientations of neighboring molecules, 
local ordering motifs for oligopyridine supramolecules on 
a surface have been successfully reproduced [10(. For 
mixtures of melamines, PTCDI, and PTCDA structural 
stability diag rams of two competing patterns have been 
reported [11 0J| using effective pair energies between 
neighboring molecules on a lattice. For the case of 
flexible organic molecules, we have recently shown that 
a multi-method approach combining MM, DFT, and 
Monte Carlo simulations is capable of providing reliable 
predictions for the emerging multiple coexisting patterns 
on surfaces fl3l4l5B. 



In this paper we provide a detailed discussion of the 
theoretical aspects of the use of interaction-site models 
for Frechet dendrons. In particular, we discuss how sa- 
lient features of the molecular structure can be identi- 
fied and used to tailor a suitable coarse-grained model 
and then be studied by means of Monte-Carlo methods. 
Furthermore, we show that a phase diversity naturally 
emerges in agreement with experimental studies and dis- 
cuss the corresponding phase stability as the packing 
fraction on the substrate or the temperature is varied. 
Then we provide a thorough discussion on the merits of 
the approach and its limitations as well as possible ex- 
tensions. 
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II. EXPERIMENTALLY OBSERVED PATTERNS 

In this Section we first summarize the experimental 
findin gs o n molecular self-organization for Frechet den- 
drons jlrl Il7| . These second generation Frechet dendrons 
are flexible supramolecules, which consist of three phenyl 
rings symmetrically disposed at the vertices of a triangle 
with two alkoxychains attached at each lateral phenyl 
ring. The arms of these chains consist at one side of the 
molecule of twelve carbon atoms, and on the other side of 
eight carbon atoms — see the structural formula in Fig. Q] 




Figure 1: Stick model of the asymmetric Frechet dendron used 
in the experiments. 

Scanning tunneling microscopy (STM) images show 
that a large variety of self-organized ordering motifs 
on top of highly oriented pyrolytic graphite (HOPG) 
emerges fl8l . Il9| . This substrate exhibits a honeycomb- 
like surface with a six-fold rotational symmetry. Seven 
general ordering motifs have been reported, five of which 
are flat-lying configurations [j^ . The tiretrack, wave, 
honeycomb, sawtooth, and jigsaw motifs are displayed in 
Fig. [21 These patterns coexist on the HOPG surface and 
phase transform over time into the thermodynamically 
stable tiretrack pattern. The schematics of the ordering 
motifs, Fig. O have been obtained employing atomistic 
modeling by molecular mechanics energy minimizations 
using the Forcite module of Materials studio 4.3 and em- 
ploying a universal force field. 



III. THE INTERACTION-SITE MODEL 

To obtain reliable predictions of the emergent pat- 
tern in the assembly of monolayers without any a pri- 
ori knowledge from the experiments, a drastic reduction 
of complexity is highly desirable. Here, we propose an 
interaction- site model which aims to provide a suitable 
description of the self-organization process. The basic 
idea is to reduce the numerous microscopic forces to a 
few representative interactions acting on groups of selec- 
ted points. The major challenge consists of appropriately 
identifying the positions of the interaction centers and 
forces relevant for the self-organization. We will demon- 
strate in this paper that the geometry of the building 
blocks, the symmetry of the substrate, and the coupling 
of both are key elements in the self-organization of mono- 
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Figure 2: Five experimentally observed patterns for Frechet 
dendrons. STM images on the left and the corresponding 
schematic ordering motifs on the right. 



layers. We will show for an experimental model based on 
Frechet dendrons, that the interaction-site model prop- 
erly captures the essence of the system, predicts the ob- 
served patterns, as well as the temperature regime at 
which they are stable. Moreover, monitoring suitable or- 
der parameters, we are able to identify in principle the 
nature and location of the melting transition. 



A. Construction of the interaction sites 

The Frechet dendron introduced above is a suitable 
and especially versatile model system to asses the validity 
of an interaction-site approach in the prediction of self- 
organized monolayers. According to the experimental 
observations the crucial features leading to the self- 
organization are the steric repulsion between the molecu- 
lar rings, the weak interactions of the carbonated chains, 
as well as the coupling of the building blocks with the 
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substrate. We model the symmetric molecular core as 
three hard spheres of radii r r located at the vertices of a 
flattened isosceles triangle — large spheres in the sketch of 
Fig. [31 The long base of the triangle is twice of its height 
I, hence all three interaction sites have equal distance I 
from the center of the base. The carbonated chains of the 
molecules are modeled by a small number of sites with 
Lennard- Jones interactions. Four neighboring CH2 units 
of the alkoxychain are coarse-grained to one Lennard- 
Jones site — small spheres in Fig. [3J The interaction sites 
are arranged in straight, rigid arms with their centers 
separated by a distance a. 

The radii of the hard spheres and the range of Lennard- 
Jones interaction, the distance of the former to the mo- 
lecular center, as well as the length and orientation of the 
arms are the parameters of the model. These geomet- 
rical features are obtained from the conformations which 
minimize molecular mechanics simulations (MM) start- 
ing from an initial configuration derived from a detailed 
analysis of the experimental STM images, Fig. [2] The 
relevant length scales of the Frechet dendrons range ap- 
proximately from 15 A of the skeleton to 50 A for the ex- 
tended molecule, with I = 6.1 A, r r = 2.6 A and a = 6 A. 



the graphite surface, and the total adsorption energy is 
about ten times larger than the total intermolecular in- 
teraction. In addition, as the size of the entire molecule, 
around 45 A, is up to ten times larger than the lattice 
constant (a few A), only the symmetry of the under- 
lying substrate plays a role in the monolayer assembly. 
In our modeling the molecules are fixed at the sites of 
a coarse-grained, fully occupied triangular lattice. The 
lattice constant a is comparable to the size of the build- 
ing blocks and the lattice exhibits the same symmetry as 
the original graphite honeycomb structure. The Frechet 
dendrons may rotate by discrete angles as rigid bodies 
around their centers and adopt one of the six preferred 
orientations of the underlying graphite. 

The interaction-site model accomplishes a significant 
reduction of degrees of freedom, setting the flexibility of 
the molecule aside. While the actual physical system 
contains hundreds of atoms per molecule able to displace 
and rotate independently, the coarse-grained interaction- 
site model consists of a rigid object with no other degrees 
of freedom than the rotation around its center. 

In addition to the conformation closest to the experi- 
mental findings, Fig. [31 we explore the Frechet dendron 
in five alternative conformations by varying the orienta- 
tions of the arms, as displayed in Fig.[H where the angles 
follow the symmetry directions of the substrate. 
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Figure 3: Interaction-site model for a Frechet dendron in 
conformation a. Large spheres account for aromatic rings, 
whereas the small spheres represent subunits in the carbon- 
ated chains. Every sphere represents four CH2 units in the 
arms of the dendrons. 

We encode the strong intermolecular steric interaction 
between the aromatic rings by hard-sphere repulsions 
between the three central spheres of different molecules, 
preventing the cores of the molecules from overlapping. 
The weak, short-ranged van-der-Waals attraction of the 
lateral chains is described by a Lennard- Jones potential: 
V(r) = 4e [{cr/r) 12 — {cr/r) 6 ] among beads in chains of 
different molecules. Thus e corresponds to the minimum 
energy which occurs at r = 2 1 ' 6 cr. 

The atomically flat graphite surface which constitutes 
the template for the pattern formation, offers six ener- 
getically equivalent orientations for the molecule. The 
molecule-substrate attraction is mainly mediated by tt- 
interactions between the phenyl-rings and the chains with 
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Figure 4: Molecular conformations showing different orient- 
ations of the lateral straight arms. The angles of the arms 
with respect to the positive x-axis are given clockwise from 
the left to the right for the five additional conformations: /3 = 
(-5tt/6, 5tt/6, 5tt/6, -tt/6), 7 = (-5tt/6, it/2, it/2, — tt/6), 
5 = (-57r/6,7r/6, 71-/2,71-/6) e = (-5tt/6, it/2, it/2, tt/6), C = 
( — 57r/6, 7r/6, 7r/6, —tt/6). 



B. Predictions of the model 

To find the regular patterns emerging in the self- 
organized monolayers we have run Monte Carlo (MC) 
simulations [2(J - [22( | considering some hundred to a few 
thousand molecules which corresponds also to the exper- 
imental situation: a single dendron covers a surface of 
approximately 4 nm 2 and the samples imaged with STM 
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occupy some hundred nm 2 . The lattice constant which 
essentially fixes the packing fraction of the monolayer has 
been varied from a — 2.8 a to a = 4.2 a covering the ex- 
perimental regime. To ensure that the stable patterns are 
assumed in the simulations, simulated annealing (23l - [25j 
has been employed. Random configurations, where mo- 
lecules assume a random direction among the six possible 
orientations, have been chosen as initial configurations. 
Starting at a given temperature the system evolves via 
Monte Carlo moves by performing 7r/3 rotations. We 
have run 500 MC sweeps [HJ to equilibrate the sample 
and measure the energy and order parameters during an- 
other 1000 Monte Carlo sweeps before lowering or raising 
the temperature. Then the temperature is lowered and 
the process is repeated. The system is cooled down to 
temperatures around 170 K, low enough for our purpose, 
given that the experiments were performed at room tem- 
perature. 
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Figure 5: Ordered motifs found by simulated annealing for 
the Frechet dendrons with different conformations within the 
interaction-site model. Molecules are represented as wedges 
matching their three aromatic rings. The carbonated arms 
are omitted for clarity. 

An abundant variety of ordered patterns arises by cool- 
ing samples of the interaction-site model in all conforma- 
tions, see Fig. Qp T hese patterns display different degrees 
of symmetries [26| which are subgroups of the p6m sym- 
metry (six-fold rotation, mirror symmetry) of the under- 
lying lattice mimicking the graphite substrate. The hon- 
eycomb, trimer and inverted honeycomb structure, see 
Fig. 03 display a p6 symmetry, i.e. the sixfold rotational 
symmetry is still present, yet the mirror symmetry of 
the substrate is broken. The sawtooth pattern exhibits 
only a p3 with threefold rotations of the unit cell. In the 
tiretrack/wave pattern rows of aligned units emerge with 
opposite directions resulting in a p2gg symmetry, i.e. a 
twofold rotational symmetry and glide- mirror axes. The 
class of p2 symmetric patterns consists of the jigsaw, zip- 
per and rhomboid ordering. Here the jigsaw unit cell is 



comprised only of two units arranged in a head-to-head 
configuration, whereas the more exotic zipper contains 
four molecules in a unit cell, the rhomboid pattern even 
six building blocks. The lowest symmetry cm consisting 
of a mirror axis and a glide mirror axis, yet no rotational 
symmetry, is realized in the crown pattern, where four 
molecules order in a head-in fashion. 

Table U lists the patterns found in the cooling process 
for the different conformations. Some patterns turn out 
to be only metastable and the ground state may change 
upon varying the packing fraction. The configuration a 
which resembles most closely the experimental configur- 
ation displays a significantly larger variety of patterns. 

To asses the stability of the emergent patterns, we have 
prepared perfectly ordered configurations at low tem- 
peratures and heated them up slowly while monitoring 
the evolution of the energy and suitable order paramet- 
ers. Thus, one can identify metastable phases which dir- 
ectly transform into a disordered state. In particular, 
one can estimate the melting temperature of the broken 
symmetry phase and determine whether the transition is 
continuous or discontinuous. 

The order parameter is defined explicitly for every pat- 
tern, based on the fraction of molecules in a given sub- 
lattice following its preferred orientation. Indeed, one 
can distinguish various sublattices in which the molecules 
point in the same direction; for instance, alternating rows 
in the tiretrack/wave pattern represent two different sub- 
lattices. Thus, for a sublattice A with Na molecules 
and preferred orientation o A , the partial order parameter 
reads 
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The factor 1/5 is introduced to ensure that ttia = in a 
disordered phase where all six orientations are equiprob- 
able. The global order parameter we monitor is just the 
average of the order parameters nil over all sublattices n a 
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The thermal fluctuations of the order parameter encode 
the linear response of the system with respect to a fic- 
titious external aligning field. Then a susceptibility can 
be defined as 
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and its behavior as a function of temperature is indicative 
of the nature of the transition. Similarly, we measure 
in the Monte Carlo simulation the average energy per 
molecule (u) and the corresponding fluctuations 
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which represents the specific heat per particle. 

The simulation results are exemplified in Fig. [5] for the 
melting transition of the sawtooth phase as a function of 
the reduced temperature fc^T/e. The average energy and 
the order parameter drop in a small temperature inter- 
val suggesting a first order melting transition. The cor- 
responding susceptibilities exhibit corresponding peaks 
close to the transition temperature, and are interpreted 
as smeared delta functions. Yet, we cannot exclude a 
continuous transition, and finite size scaling would be 
needed to clarify the order of the transition. In this work 
we focus on the pattern diversity and have not pursued 
this issue any further. 
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Table I: Emerging ordered phases by cooling down samples 
with N — 576 Frechet dendrons for the chain conformations 
a, . . . , f . The phases in bold font are stable under heating, 
while the other phases would be unstable in an experiment at 
room temperature. 




Figure 6: Melting transition for the sawtooth phase. Top: 
Energy and heat capacity per molecule as a function of the 
dimensionless temperature. Botom: Order parameter and 
suscecptibility. The simulation was performed for N = 576 
Frechet dendrons in the /3 conformation on a lattice with 
a/a = 3.8. 



C. Comparison to experiments 

The interaction-site model reproduces many of the fea- 
tures of the experimental model system consisting of 
Frechet dendrons on graphite surfaces. First, in both 
systems an ample diversity of patterns emerges suggest- 
ing that the basic building blocks are correctly trans- 
ferred from the Molecular Mechanics calculations to the 
interaction-site model. Furthermore, the patterns found 
in the theoretical model closely resemble the experi- 
mental ones both in its global and local ordering motifs, 
i.e. the structure displays the same wallpaper group (26| 
and a similar arrangement of the molecules in a unit cell. 

In the highly symmetric trimer honeycomb structure 
three molecules align in a subunit facing each other. The 
resulting trimer exhibits a threefold rotational symmetry 
without additional mirror symmetry implying that the 
trimer is chiral. These trimers order in sixfold symmetric 



arrangement where a single molecule of unspecified ori- 
entation resides at its center, see Fig.Ek). In the experi- 
mental system precisely the same structure has been ob- 
served where the center of the hexagons appear as blurred 
regions in the STM measurements (l5| . Thus it appears 
that these unpaired molecules are free to change their 
orientation rapidly, much faster than the time resolution 
of the STM (23. 

The tiretrack/ wave pattern is characterized by altern- 
ating columns of molecules pointing in the same direc- 
tion. Then in the corresponding rows of this pattern 
every second molecule possesses the same orientation. 
Depending on the conformation of the arms of the mo- 
lecule, the pattern is formed by strong intracolumn inter- 
actions with moderately coupled columns, or by strongly 
linked rows that then arrange in a parallel fashion. Sim- 
ilar to the theoretical model, the experimental system 
displays a tiretrack and a wave pattern where molecules 
arrange in rows or columns, respectively, see Fig.^J. Ex- 
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perimentally the coupling in the rows and columns is sig- 
nificantly different resulting into two clearly distinguish- 
able phases, although wave and tiretrack belong to the 
same wallpaper group. For the simulations on a coarse- 
grained lattice this distinction can no longer be made and 
both phases merge into a single tiretrack/wave pattern. 

Our simulations also reproduce the sawtooth pattern, 
Fig. [5]i, which has been observed in STM measurements 
on Frechet dendrons. Here three molecules form a com- 
posite which constitutes the repeat unit on a triangular 
lattice. In contrast the experiments report a lower wall- 
paper group, p2 rather than p3, since here two trimers 
arrange in a opposite orientation to form a unit cell, com- 
pare Fig. [2J 

For the case of the jigsaw pattern, Fig. [Sf, we find 
almost perfect agreement with the STM images, both 
with respect to the local ordering as well as with respect 
to the wallpaper group. 

The remaining patterns we have generated in the 
interaction-site model have not been found for Frechet 
dendrons. The first group of patterns (Fig. [5k,-c) are 
all highly symmetric and very similar in their respective 
local motifs and therefore these structures may be very 
sensitive to the conformation and the details of the mo- 
lecule. Thus it appears promising to modify the chemical 
structure of the Frechet dendron only slightly to realize 
also the inverted honeycomb and the honeycomb struc- 
ture. The second group of simulated patterns without 
corresponding experimental result consists of large com- 
plex unit cell with low symmetry (Fig. Eg-i). Therefore 
even if some of these patterns constitute the ground state 
of the system it is likely that they are not realized exper- 
imentally due to kinetic barriers. 
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Figure 7: Zero temperature energies for conformation 7 as a 
function of the lattice constant. The ground state correspond- 
ing to the lowest energy depends sensitively on the packing 
fraction demonstrating that the competing ordering motifs 
display similar energies. 

The energies of the various patterns are rather sim- 
ilar in agreement with the experimental observation of 



coexisting metastable phases. After annealing in the ex- 
periments or sufficiently many Monte Carlo steps in the 
simulation the patterns are expected to transform to the 
stable pattern. For temperatures relevant in the experi- 
ment these patterns are almost always identified with the 
ground state of the system. Which of the suggested pat- 
terns actually represents the ground state depends critic- 
ally on the density. We exemplify the subtle interplay of 
the effects of packing and ordering for conformation 7 in 
Fig. [7]which displays one of the most complex behaviors. 
For high packing fractions (small lattice constant) first 
the honeycomb pattern and then the sawtooth pattern 
corresponds to the ground state. For moderate densities 
the tiretrack/wave exhibits the lowest energy and exper- 
imentally it also is found as the most stable pattern. 

We have shown that interaction-site models are a 
powerful approach to model pattern formation of com- 
plex supramolecules on a substrate. In particular, the 
pattern diversity is correctly transferred from the mo- 
lecular mechanics minimizations and also details of the 
patterns such as symmetries and local ordering motifs 
are surprisingly well reproduced. Let us discuss limita- 
tions and possible improvements of the current model. 
First, the interaction-site model so far accounts only for 
the dominant van der Waals interactions, yet on closer 
inspection electrostatic interactions could be also import- 
ant for the selection of a particular motif. From the mo- 
lecular mechanics calculations one infers that these forces 
are responsible for up to 30% of the intermolecular inter- 
action. It would be desirable to model also these elec- 
trostatic interactions and test if more accurate results 
can be achieved. Similarly, in the case of highly direc- 
tional forces, such as hydrogen bonds, one should expli- 
citly incorporate these specific interactions in the inter- 
action site model. Second, the model does not account 
for the effects of the solvent. Nevertheless some solvents 
interact via hydrogen bonds with the Frechet dendrons 
and may favor conformations of the Frechet dendrons 
where oxygen atoms are exposed to the solvent. Since 
the conformations are crucial for the corresponding pat- 
terns one should also account for the solvent effect in the 
interaction-site model. For example, one could attribute 
an additional conformational energy for each molecule 
and then compare the total energy, consisting of both 
conformational as well configurational energy, and thus 
determine the stable pattern and stable conformation in 
a specific solvent. Furthermore the lattices considered 
so far are fully occupied, which excludes complex motifs 
with host structures. In the interaction-site model one 
could account for this possibility by either a fixed partial 
coverage density (canonical ensemble) or a fixed chemical 
potential (grand canonical ensemble) and include hop- 
ping processes in addition to the rotational Monte Carlo 
moves. Last, we have reduced the molecule-substrate in- 
teraction to the lattice symmetry and spacing. The ma- 
jor justification for this simplification is that the supra- 
molecules are significantly larger than a unit cell of the 
substrate and correspondingly experience only an aver- 
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aged adsorption energy with preferential directions. Yet, 
for smaller molecules, larger substrate lattice constants, 
or specific interactions one should include explicit space- 
resolved molecule-substrate interactions. 



IV. CONCLUSIONS AND OUTLOOK 

In this work, we have developed an interaction-site 
model which correctly predicts the ordered motifs of as- 
sembled monolayers. By reducing the degrees of free- 
dom and considering the building-blocks as rigid bodies 
with a reduced set of interacting points, we have demon- 
strated that the self-organization relies on very general 
features of the system considered: the coupling with 
the substrate, the geometry of the building blocks, and 
their weak interactions. These are universal principles 
in self-organizing systems which do not depend on the 
specific nature of the building blocks and the underlying 
substrate — our method also works for substrates exhibit- 
ing different geometries. Therefore, the predictive power 
provided by our model may guide the synthesis of suitable 
building blocks to engineer arbitrary patterns for specific 
goals. The versatility offered to construct the building 
blocks makes our model especially suitable to explore a 
wide range of geometries. 

In addition, we have shown that by combining sim- 
ulated annealing with a subsequent slow heating of the 
system, one can not only predict the emergent patterns, 
but also their stability upon heating and the nature of the 
transition into a disordered phase. We have found that a 
broad variety of long-range ordered phases are stable for 
various conformations of the building blocks and density 
regimes which may indeed coexist, as it has been ob- 
served in the experiments. The melting temperatures of 
the ordered motifs range from approximately 500 K to 
1500 K, much higher than the room temperature where 
the experiments were performed. 



However, the model still lacks an active determination 
of the intramolecular conformation. Here we have con- 
sidered the building blocks as rigid bodies with the con- 
formations observed in the molecular mechanics minim- 
izations. A natural generalization would be to make the 
building blocks flexible and explore the interplay between 
intra and inter-molecular ordering. This extension im- 
plies including at least four new degrees of freedom per 
molecule, the orientations of the arms, which makes the 
Monte Carlo simulations computationally very expens- 
ive. An interesting alternative to speed up simulations 
would be to rely on genetic algorithms (28l . |29| . They are 
suitable methods to compute the ground states of two- 
dimensional systems by minimizing the energy of single 
unit cells [3(1 [H|. The conformations and motifs res- 
ulting from this minimization can serve as input config- 
urations for the Monte Carlo simulations to assess their 
stability. In this way, the interaction-site model is less 
dependent on the details of the external input. In ad- 
dition, when employing genetic algorithms one can relax 
the constraint of the substrate. Instead considering the 
building blocks to be attached to the sites of a lattice, 
one can mimic the role of the substrate through arbit- 
rary potentials. This opens a way to investigate systems 
whose substrates display more complex symmetries and 
interactions with the building-blocks. 
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